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Abstract 

Within some reasonable approximations we calculate deformation of 
an elastic bar, falling on a source of central gravitational field. We consider 
both elastic deformations and plastic flow together with destruction of the 
bar. Concrete calculations for a number of materials are presented. 
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1 Introduction 

The study of the extended deformable bodies motion in general relativity (GR) 
is an important and technically challenging problem with fundamental as well 
as practical applications. In order to be consistent with principles of GR the so- 
lution requires a thorough review of all the notions which are commonly used in 
nonrelativistic continuous media physics. General principles lying at the foun- 
dation of deformable bodies physics were actively debated when the relativity 
theory first appeared [TJ [2J [3]. Contemporary understanding of the problem 
is based upon seminal works [H [5j [6] . Different aspects of the elasticity the- 
ory in GR were further developed in works [7] (gauge formulation), [5] and 
[HI EH E] (general theorems of existence and uniqueness) . The most important 
and widely accepted topics in the field are: material (frame independent) de- 
scription of elastic medium, 1+3- decomposition of space-time and Lagrangian 
formulation of gravitoelastic equations. 

Note that more detailed treatment and application of these topics gave rise 
to slightly different versions of the general approach, depending on author's 
conceptions; some of these conceptions are outside the scope of elasticity in GR 
[T21 IT51 [18] . Despite that these slightly different approaches to the problem 
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reproduce the right nonrelativistic limit, the predictions regarding experiments 
with extended elastic bodies based on them can diverge drastically. That applies 
to the behavior of the bodies in strong gravitational fields (horizon of black hole) 
and even in weak fields without Newtonian limit (gravitational waves) . So, the 
study of the behavior of deformable bodies in such exotic situations (from the 
point of view of Newtonian gravity) could serve for improvement and further 
development of elasticity theory in GR. At present both gravitational waves 
physics and black hole physics are theoretically well founded and actively devel- 
oping topics. However, while the former has many intersections with elasticity 
in GR (first of all bar detectors and their interaction with wave perturbation 
in the space-time metric should be mentioned) , incorporation of elasticity ideas 
into the black holes physics is much more modest. The majority of attempts 
in this or related fields dealt with the self-consistent problem of equilibrium of 
astronomical bodies that possess elastic properties (among the last works on the 
subject is already cited [13]) and with tidal phenomena in celestial mechanics 
and astrophysics [H j [T5 | [To ] IT?] . 

The present paper can be treated as a small step towards understanding what 
will be the influence of a strong gravitational field on test deformable bodies. 
Note that this problem is not the subject of traditional elasticity GR since 
strong gravitational field causes tidal stresses which could lead to plastic flow and 
mechanical destruction of any real body. We are going to illustrate our approach 
by a relatively simple example of a one-dimensional continuous medium that is 
thin elastic bar, falling towards the source of central gravitational field. We'll 
show that our treatment of the problem encounters two difficulties which can 
be analyzed independently. The first difficulty originates from nonlinearity of 
the equations of motion and has nothing to do with relativistic aspects of the 
problem whatsoever. The second difficulty arises from dealing with relativistic 
properties of both motion and space-time that have to taken into account in 
order to get a correct quantitative picture of deformations. 

The aim of this part of the paper is to solve the first kind of difficulties within 
Newtonian gravity and to lay a foundation for more precise and consistent con- 
sideration in the context of GR. We formulate an approximate "quasistatic 
fall" approach which allows to solve nonlinear equations by the perturbations 
method, with smallness parameter L/r, where L is the undeformed bar's length, 
r is instant distance from the bar to the source of gravitational field. The law 
of motion of point-like particles will play the role of " null approximation" . This 
nonrelativistic approach allows us to calculate the deformational characteristics 
of bars up their destruction. In section [5] we derive the exact equation of mo- 
tion (within Newtonian theory) that describes the bar's behaviour including the 
deformation. In section we discuss the exact formulation of the quasistatic 
fall condition. We find out that it is valid for bars made of solid materials. In 
section UJ we derive an approximate equation of motion in first order on L/r 
and find solution to the derived equation that satisfies all necessary initial and 
boundary conditions. Our approximate solution for elasticity reproduces the 
exact static solution (see [19]) after the transversal averaging and coarsening 
procedures. Quantitative and qualitative deformational characteristics of the 
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bar are calculated in section In section [5] we derive slightly simpler formulae 
for calculation of definite deformational characteristics: the time of the begin- 
ning of the plastic flow and the time of destruction as well as the corresponding 
positions of the bar. Then, using the derived formulae, we find some "universal" 
relations independent of the bar's initial state and present results acquired for 
some specified materials. In the next part of the paper there will be presented 
relativistic analysis of the problem considering a deformable body falling into 
Schwarzschild black hole. 

2 Statement of the problem 

Let us consider radial fall of a thin probe bar from an initial position distanced 
from a massive source of Newtonian central gravitational field. While in motion 
the bar stays oriented along the line of force (in a radial direction). Let us 
put the origin O of an one-dimensional coordinate system at the force center 
and direct the coordinate axe OX along the bar's motion radial line. Then 
kinematics of the bar can be described by a displacement field x(t, £), defined 
on the points £ of a remote unstrained bar (Lagrangian picture j]], and depending 
on time. The initial and boundary conditions are as follows: 



i(0,0 =ro + £, £e [-L/2-L/2]; x(0, £) = -vq = const; ct| ?=0 , l = 0. (1) 



Here ro is the initial coordinate for the bar's center of mass, L is the unstrained 
bar's length, vq is absolute value of initial velocity, a — a(t,£) — a xx is 1- 
dimensional "stress tensor" (later to be referred as "stress") that depends on £ 
and t. Here and below (if not specified otherwise) the dot denotes differentiation 
with respect to time coordinate t (or r — see below) , the accent denotes differ- 
entiation with respect to space Lagrangian coordinate £. To be more precise, 
we define a as the following average value: 



where integral is being taking over transversal section of the bar, S means its 
area and cr xx is true longitudinal stress tensor component, that depends on 
the points of the transversal section. 

Lets us consider small element dx of the bar at some fixed moment of time 
t. Its ends have material coordinates £ and £ + d£. Inertial force, related to the 

1 So, we absolutely ignore transversal elasticity of the bar. Such simplifying is correct, 
since ratio of absolute transversal Ar and longitudinal Ai deformations satisfies the following 
nonequality: 



due to thinness of the bar (R is its characteristic transversal size) and smallness of Poisson's 
coefficient for majority of solid materials. Exact results of [19] support this assumption: even 
for solid bar with R/L = 1/4 ratio Ar max /AZ ~ 1/10 near the Earth. 
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element is expressed by the formula: 



—dm'x = —pSxdx, 

where p = p(t, £) is the local volume mass density of the bar. Gravitational 
force, acting on the same element, has the form: 

GMdm GMpSdx 



where M is mass of origin, G is the Newtonian constant. Finally, elastic force, 
acting on the element looks like this: 



d(aS) 
dx 



dx. 



After we've added the three forces together and equated the resultant to zero, 
the following equation of motion can be obtained: 



1 d{aS) 
pS dx 



GM 



(2) 



In order to exclude stress it is necessary to set some defining relation [5D] for 
the material of the bar. Said relation normally connects stress a with relative 
deformation e = x' — 1. In this paper we'll use piecewise linear defining relation, 
which is commonly applied in the majority of the model problems dealing with 
elasticity theory and strength of materials [21j . FigQ]shows a simplified effective 
diagram "strain-stress" illustrating this defining relation. In real laboratory 



vS/So 



£1 



Figure 1: Simplified stress diagram of bar's material, expressing results of laboratory 
experiments with real bar. 



experiments with bars the controlled external force is usually related to the 
initial area So of the transversal section before stress starts to build up. So, such 
experiments allow only to determine dependency of aS/So on e, which includes 
variation of the transversal section area under deformation. The inclined part of 
the diagram describes elastic deformations (Hooke 's law), for which aS/ Sq = Ee, 
where E = O0A0 is Young's modulus for the material of the bar, eo is limit of 
elasticity, Co is stress of plastic flow. The latter is described by the horizontal 
part of the diagram, where strain increases under practically constant stress. 
Parameter ei characterizes ultimate relative deformation causing the destruction 
of the material. 
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Using the chosen defining relation, the term linked to stress in © can be 
transformed (under elastic deformations) as follows: 



So 



d{Ee) 



So 



(E(x'-l)) 
dx/dS, 



dx 



where in the differentiation with respect to x we've turned to the Lagrange 
variable £. Assuming that the bar is homogeneous, we have E = const and 
po = /j(0,£) = const. Using of the relation: 



that expresses the law of mass conservation as applied to any element of the bar 
related to fixed material particles, equation ([2]) can be rewritten in the form: 



where □ = d 2 — <9| is the two-dimensional D'Alembert operator, r = cot, Co = 
W E I po is sound velocity in the material of the bar, r s — 2GM/cq is sound 
analog of gravitational radius, which we' 11 call the sound radiu& 

Equation ([3|) is nonlinear wave equation, for which standard methods of 
mathematical physics are inapplicable. It easy to find a family of particular 
soliton-like solutions, which is characterized by special dependency on £ and 
r: x = x(Q, £ = ar + where a and b are arbitrary constants, satisfying 
the condition: \a\ ^ Substitution to ([3]) leads to an ordinary differential 
equation: 



with k = r s /2(a 2 — b 2 ), where the accent denotes differentiation with respect to 
£. However, its general solution: 



--^(Ax + 2k)x+ — -ln(k + Ax+ v / Ax(Ax + 2k)) = (- Co, (5) 



where A and Co are integration constants, doesn't satisfy conditions ([T]), no 
matter what constants or parameters are being chosen. In general, our study 
will not involve those stressed states of the bar that are associated with running 
waves. We'll restrict ourselves to seeking solutions to the equation © that 
describe a smooth (quasi-static) variation of the bar's stressed state from null 
strain (far from the gravitational center) up to the point of destruction. Note, 
that for those parts of the bar, where elasticity limit is exceeded, the equation 
((3]), according to the stress diagram in Fig[IJ will acquire a simpler form as that 
of an ordinary differential equation similar to ^ . 

2 In a difference with gravitational one, it depends on elastic properties of bodies, falling 
on the gravitational center. 



PoSq = p(t,£)x'S, 




x- 



(4) 



A- 



5 



3 Condition of quasi-static fall 



On an obvious assumption that the law of motion for material points roughly 
applies to the motion of extended small-sized bodies, an approximate solution 
to the equation ([3]) can be constructed. To be more precise, we state that: 

lim max \x(t) — xr\(t)\ = 0, (6) 

L^0t G [0;T] 

where L is characteristic size of the extended body, x(t) is exact law of motion 
for the body's center of mass, calculated for time interval [0;T] with regard 
to its extended structure, xo(t) is law of motion for point-like body with the 
same integral characteristic (mass, charge etc.), calculated for the same time 
interval. The relation © is, in fact, one of the basic statements of both classical 
mechanics and relativity theory in implicit form. (Precise formulations and 
important rigorous theorems see in [TO]). 

This assumption allows us to seek solution to the equation as follows: 

x(t,0=x (t)+£+x(t,0, (7) 

where xq (t) + £ describes the motions of the rigid bar (tidal forces to be dis- 
regarded), which center of mass initially is in the position xo(0) = ro, while 
correction x( T i£) describes strain of the bar, subjected to in some set mode 
of motion Xq(t). After we've calculated the strain in a set mode of motion, if 
necessary, we can define the mode of fall more accurately and recalculate the 
strain etc. 

In order to establish the limits of validity of our approach, let us estimate 
critical distance r c , over which the bar is destroyed. Upon equating the char- 
acteristic tidal stress to ultimate stress, which is for majority of solid materials 
three orders less then their Young's modulus, our estimation appears as follows: 

r c ~10(r g L 2 ) 1/3 (c/c ) 2 / 3 , 

where r g = 2GM/c 2 is the gravitational radius of the field's source. Assum- 
ing that c/cq ~ 10 5 , we obtain a more convenient formula for calculating the 
estimations: 

r c ~ lOVgi 2 ) 1 / 3 . (8) 

It follows from ([8]) that for gravitational centers whose mass is M < 10 3 M Q and 
bars of the lengths L ~ 10°m r c > r g (for the Sun r g ~ 1km and r c ~ 100km). 

The validity condition of the approach discussed above is the "quasi-static 
fall" quality: time AT for characteristic variation of the tidal force must exceed 
by far the time L/cq for the longitudinal sound waves propagation along the 
bar. The time AT can be roughly estimated as xq/coXq- As the expression 
implies, the time will be minimal at the moment of destruction. The energy 
conservation law has the following dimensionless form: 

^ = - + 4-- (9) 

x eg r 
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The first term on the right at the moment of destruction has the order 10 7 
(when L ~ 10°m, M ~ M©), while the second and third under nonrelativistic 
initial velocities and large enough (compared to the gravitational radius of the 
source) initial distances have considerably less order. Thus we can disregard 
them at the moment of destruction and applying ||5J) we get: 

AT min ~10 3 / 2 ->-. (10) 
Co c 

So, for bars made of solid materials (metals, steels and alloys, for which u\ ~ 
1Q~ 3 E and Co ~ 10~ 5 c) the condition of the fall's quasi-static quality is met 
automatically by application of its law of motion. For softer materials the con- 
dition can be violated near the point of destruction. That implies the necessity 
of studying the shock waves and sound retarding impact on the problem. 



4 Approximate solution 

After we have decomposed the gravitational force in the row (with respect to 
L/xq) near the rigid bar's center of mass (its instant position) xq(t), the right- 
hand side of the equation ([3]) can be rewritten in the following form: 

r s r s r a £ ( L 



2{x + 2 ' x 3 ' " \x 0/ > ' ' 

After substituting it into (J3J) and taking ((7|) into account, we obtain the following 
approximate equation of motion: 

x o + X-X =-oZ5 + -Z3-- 

The first terms on the left and on the right are cancelled as the definition of xq (t) 
implicates. The second term on the left is responsible for back reaction shown 
by motion on the deformation of the bar. The term is of the order o(L 2 /xq) that 
can be proven later when the approximate solution for \ is obtained, meaning 
that we can disregard it for the time being. Then the equation that define 
strains, takes the form: 

x 

Its general solution is: 

X = -^C + C 1 (r)e + C 2 (r), 



G.r, 



where C%, Ci are arbitrary functional constants of integration. Solution x(t, £), 
that satisfies both initial (with accuracy £ 3 / r o) an d boundary conditions ((T|), 
takes the form: 



-^l(^-y): (!->> 
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5 Elastic and plastic deformations pictures of 
the falling bar 

Upon defining the displacements field u(r, £) by the formula: x — xq + £ + 
the solution (JTSJ) results in expression: 



(13) 



Figf5] shows the displacement diagram at some arbitrary (but preceding the 
plastic flow) moment of time. 
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Figure 2: Instant displacements diagram related to the center of the bar. it m ax = 
r s L 3 /24:xl 

The diagram during the fall is stretched in vertical direction. The instant 
effective stress diagram, expressed by the formula: 



Ceil 



aS r s E ( L 



(14) 



5 ZXq 

is shown in Fig|3] As the previous diagram, this one shows tendency to vertical 

'iff 




-L/2 L/2 t 

Figure 3: Instant effective stress diagram before plastic flow, cr max = Er B L 2 /8xfy 
stretching of the same time law. At the moment r* , defined by relation: 



stress reaches the first stress limit (To at the middle point of the bar and this 
marks the nucleus of the plastic phase. The size 2£b of the part of the bar 
engaged in the plastic phase increases in time symmetrically in both directions 
from the center of the bar according to the expression: 



(16) 




This dependency is shown in Fig. [¥] 

'^elastic phase 




Figure 4: Boundary of elastic and plastic phases in falling bar as function of time. 

When the front of the plastic phase part approaches the ends of the bar, 
its velocity d^b/dt asymptotically tends to zero. Any element of the bar in the 
plastic phase (i.e. situated inside the interval [— £t>(T), £t>( T )] a t the moment r) 
moves free in gravitational field, as the stress diagram gives a' = for any point 
of plastic phase. 

Since the appearence of the plastic phase, matter elements situated in the 
vicinity of the middle point of the bar, continue to suffer the stretch strain. 
However, they doesn't stretch the same way the parts of elastic phase do (i.e. 
in accordance to the solution (fT2|) ). but as particles of noninteracting "dust 
packet". It is in the central point £ = the relative plastic deformation will 
reach its maximum t\ at the moment r**, causing the bar to break in two halves. 

In order to calculate the relative deformation at the middle point of the bar 
at an arbitrary moment of the plastic phase, it is worth noting that the element 
d£ of the bar containing middle point, will by the beginning of the plastic phase 
acquire the length ds = (cq + 1) d£. Its further stretching can be described by 
the equation: 

(e + 1) d(. 

Here xq(t, s) means a family of solutions, that describe free particles with coordi- 
nates xq(t*) + s at the moment t — t* and velocity xq(t*). Relative deformation 
with respect to initial configuration is given by the expression: 

, _ds' dx (T,s) 

e ~ d£ ~ ds 

Then equation for r** takes the form: 

dx (T**,s) 
ds 

9 



ds' 



dx (r, s) 



ds 



ds 



dx (T,s) 



(eo + 1) - 1. 



s=0 



61 + 1 

e + l' 



(17) 



Equations (fTU)) - (fTT)) together with the solution ([5]) for free point-like particles 
give a complete picture of deformations of the bar up to its destruction. 



6 Design formulae and calculations 

Before we start deriving design formulae for calculations we want to stress that 
any more or less distinctive relative deformations of bar (say, about 0.01%) 
appears only at the latest times, when the bar moves considerably close to the 
gravitational center. So, for the source of mass about of the Sun's and the bar 
about lm long we can expect by (fT4|) the deformation e = 10~ 4 only at distance 
of about f ~ 100km from the gravitational center. The general formula is: 



It gives us practically the same result in the majority of the other reasonable 
situations. So, the largest and most important part of the bar's deformational 
history takes place under x < r. It implies, that instead of the general solution 
(O we can use its simpler asymptotic form under small x, where the law of 
motion is parabolic: 

^) (ro-r) 2 /^^j (to~i) 2/3 , (18) 

where to = r o/co means total time of particle's fall. This is an unique value, 
which depends on the initial state of the falling object. It can be calculated by 
means of the exact formula ([5]) : 



where Eq = Eq/tti is initial specific energy of the falling object, is initial 
distance from the gravitational center and the function 

4>(z) = y/z+l - — p ln(l + 2z + 2y/z(l + z)). 
2y z 

For r ~ 1 a.u., M ~ M Q and v = 50km/s the formula (fill)) gives us to ~ 25 
days. The exact validity criterion for the parabolic law of motion near the source 
is as follows: 

^0 „ fr s \ 2/3 _ l/3 



It is satisfied when the initial velocities are nonrelativistic and initial distances 
exceed f by far. 
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By using (01 and general formulae (flU|) - (fl~T)) of the previous section, it is 
easy to derive the following convenient design expressions, that describe all 
characteristic events of the bar's deformational history: 



to 



3x/27 



(20) 



Derivation of (fT7|) is slightly more complicated. Decomposition of |H) in the 
vicinity of the parabolic law has the following form: 



C-t = 



dx 



2x 3/2 Ax 5/2 



y/A + r s /x 



5r ; 



3/2 



(21) 



Here C and A are integration constants determined by initial conditions at the 
moment of the appearance of the plastic phase: 



x(t*) = X q(t*) + s; x(t*)=± (t*). 



(22) 



Upon substituting (|2"T]) in (|2"2"|) and solving equations with respect to C and A 
(in the first order with respect to s) we obtain the following implicit equation 
that defines the law of motion in s- vicinity of the center £ = after the moment 
of the appearance of the plastic phase: 



2x 3 / 2 



ST 



5/2 



T - T. 



(23) 



Here Xq = Xq(t*). After differentiating of the expression with respect to s under 
fixed r and assuming s = 0, we obtain: 



dx 



s=0 



5xf 



Substituting in this expression t(t) = tq(t), we finally get: 



dx 
ds 



4 
5 ! 



-4i- 

5 x 



(24) 



Destruction conditions (fTTj) and (|24|) allow us to come up with the following 
formula for the time of destruction: 



to 



Lrf 



where r\ is root of equation: 

4 4 5( £l + 1) 

T) = 



(25) 



T) 



eo + 1 



The dependency of all obtained formulae (|20| - (|25p on initial conditions is 
based on the parameter tq only. Some of consequences of these formulae are 
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perfectly independent of the initial conditions. For example, full time of plastic 
phase existence: 

At EE T** - T* = ^— (1 - 7f) (26) 

and the size of the plastic phase at the moment of destruction: 

2£ b (r**) = 2€ bmax = L^/\-rf 

depend only on the elastic constants for the material of the bar. Also we have 
the following universal ratio of coordinates x** = x(t**) and x* = x(t*): 

x**=r) 2 x*. (27) 

Table 1 shows results of calculations made for the bars of fixed length and 
of different solid materials. 



Material 


£0,10""% 


Cl, % 


x*/r g 


V 


x /r g 


2£b max/^ 


Ar/(r -T*) 


Aluminium 


3.1 


45 


53 


0.56 


17 


0.984 


0.82 


Iron 


8.5 


50 


40 


0.54 


12 


0.987 


0.84 


Gold 


5.0 


40 


86 


0.58 


29 


0.981 


0.80 


Copper 


5.8 


60 


54 


0.50 


13.5 


0.992 


0.88 


Platinum 


4.4 


45 


74 


0.56 


23 


0.984 


0.82 


Lead 


3.1 


50 


142 


0.54 


41 


0.987 


0.84 


Silver 


3.2 


45 


82 


0.56 


26 


0.984 


0.82 


Titanium 


9.1 


70 


39 


0.47 


9 


0.996 


0.90 


Steel 15GS 


17.5 


18 


31 


0.78 


19 


0.880 


0.53 


(small strength) 
















Steel N18K9M5T 


97 


8 


17 


0.90 


14 


0.684 


0.27 


(of high strength) 
















Steel 30H13 


75 


6 


19 


0.92 


16 


0.627 


0.22 


(martensite) 
















Alumalloy 


49 


10.5 


31 


0.86 


23 


0.772 


0.36 


(Al-Cu-Mg)D19T 
















Titanalloy 


108 


4 


20 


0.95 


18 


0.514 


0.14 


(Ti-Al-V-Cr)TS6 

















Table 1. Results of approximate calculations by formulae (|20[1 and (|25[) for central 
body with M = 2 • 10 30 kg and bars of length L = lm made of different solid materials. 
Limit of elasticity eo calculated by the formula: eo = (T0.2/E, where 00.2 is (effective) 
stress which causes residual relative deformation 0.2%. Ultimate deformation ei = 
eo + 5, where 8 is relative residual stretching after the destruction of the bar. Steels 
labels corresponds to Russian standards [22] . 

From (fl"5|) and (|2T|) it follows that immediately after destruction the plastic 
phase disappears in both halves under the condition 77 > 1/2 1 / 3 w 0.794. Fifth 
column of the Table 1 shows that for all listed materials except high alloy steels, 
the plastic phase will reappear in future generations of the bar's fragments. 
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Inside of the bars made of high alloy steel, plastic phase appears, increases, 
disappears and after the fission reappears again in fragments some time later. 
According to the first formula of (|2D)l . bars with lengths of about 10°cm or 
less begins to flow near the gravitational radius of the central body, which one 
should consider already a black hole. More precise and consequent calculations 
in the context of GR will be presented in next part of the paper. 



7 Conclusion 

We have considered intermediate nonrelativistic problem concerning the fall of 
an elastic bar on massive source of central gravitational field. Within the quasi- 
static fall approach formulae (|20|) and (|25|) we've been able to acquire reasonable 
and in some sense exhaustive description of the bar's deformational history up 
to its destruction. Note that our simple estimations have shown that for bars 
made of solid materials like solid metals, steels and alloys the quasistatic fall 
condition is satisfied automatically. Assumptions used for derivation of p0[) 
from ([9]) concern general nonrelativistic approximations of the problem. 

Results of calculations, presented in Table 1 show that bars with length 
about lm falling on source with mass about Mq are being destroyed far from 
the gravitational radius (closest (of the materials listed in Table 1) to it comes 
the titanium bar — about 10r g ). Note that the time At of the most important 
part of the bar's deformational history is about 1ms in absolute units. After 
the first destruction each fragment of the bar gets divided again etc. However, 
quantitative treatment of the behaviour of the generations of the bar's frag- 
ments becomes more complicated, because of non-elastic hysteresis and strain 
hardening phenomena, that take place in deformed fragments. Taking average 
values fj sa 0.67 and x* ~ 53r g from the Table 1, and using the universal formula 
([27) we obtain: 

N ~ log^ ' 



— a rough estimation of number of the consequent fragments continuing dividing 
until they reach the gravitational radius. Near the gravitational radius there 
will be required a cardinal reconsideration of our formulae within the context 
of GR. 

It is interesting, that the bar, falling in central gravitational field is, in some 
sense, a perfect experimental device for obtaining the exact stress diagrams. In 
standard laboratory experiments with bars initial homogeneous stressed state 
of the bar near ultimate stress becomes complex nonhomogeneous one with 
spontaneously appearing flowing neck. Unlike these, inside falling bar we have 
the tidal gravitational stretching that reaches its maximum at the center of the 
bar at any moment of time. Due to the volume rather than surface character 
of tidal forces, the stress diagram obtained by experiments with falling bars, 
would be independent of the end's stressing and Saint- Venant principle. 
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